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Abstract 

Wc use the available two diflForont self-consistent formulations of quasiparticle models and extend 
their applications for the description of quark gluon plasma (QGP) at non-vanishing baryon 
chemical potentials. The thermodynamical quantities calculated from these models are compared 
with the values obtained from lattice simulations and a good agreement between theoretical 
calculations and lattice QCD data suggests that the values of the parameters used in the paper 
are consistent. A new equation of state (EOS) for a gas of extended baryons and pointlikc 
mesons is presented here which incorporates the repulsive hard-core interactions arising due 
to geometrical size of baryons. A first order deconfining phase transition is constructed using 
Gibb's equilibrium criteria between the hadron gas EOS and quasiparticle model EOS for the 
weakly interacting quark matter. This leads to an interesting finding that the phase transition 
line ends at a critical end point beyond which a crossover region exists in the phase diagram. 
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1. Introduction 

The exploration of the QCD phase diagram, both experimentally and theoretically, 
is one of the main thrusts of the present day research in strong interaction physics. 
The existence and the position of the critical end point (CEP) are the key questions 
in mapping the QCD phase diagram. Lattice gauge calculations first revealed that the 
transition between hadron gas (HG) and quark gluon plasma (QGP) phase at /Ub = 
and large T is a cross-over transition and there were further indications that the cross- 
over transition turned into a first order chiral phase transition for nonvanishing and 
finite values of /x_b [!]• Although, the existence of a CEP in conjectured phase diagram 
of strongly interacting matter was proposed a long time ago [2] , the absence of the CEP 
was also noticed in some recent lattice calculations [3]. Thus, the existence and location 
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of the CEP is still a matter of debate. Therefore, it is worthwhile to investigate the 
existence and/or non-existence of CEP in various phenomenalogical models. Moreover, 
if the CEP exists then it is also important to find the precise location of the CEP 
and to determine its properties. In this contribution, we consider QGP as a system of 
quasiparticles, which arc quarks and gluons possessing mcdium-dcpcndcnt masses arising 
due to vacuum interactions [4-7]. Furthermore, we obtain a new EOS for the hot and 
dense HG as formulated in one of our publications [8] , and construct the phase boundary 
by equating the QGP pressure with that of HG pressure using Gibbs' equilibrium criteria. 
A more detailed description can be found in Ref. [9]. 



2. EOS for a Hadron gas 



Recently we have proposed a thermodynamically consistent excluded volume model 

for hot and dense hadron gas (HG). In this model, the pressure for the HG with full 
quantum statistics and after incorporating excluded volume correction can be written as 
[10]: 

i i 

where R is the fractional occupied volume by the baryon and Aj is the fugacity of the ith 
baryon. Furthermore, in Eq. (1), the first term represents the pressure due to all types 
of baryons where excluded volume correction is incorporated and the second term gives 
the total pressure due to all mesons in HG having a pointlike size. In this calculation, 
we have taken an equal volume V° = for each type of baryon with an identical 
hard-core radius r = 0.8/m. We have taken all baryons and mesons and their resonances 
having masses upto 2GeV/c? in our calculation for HG pressure. Wc have also used the 
condition of strangeness conservation. We want to stress here that the form of this model 
used under Boltzmann approximation has been found to describe [11-12] very well the 
observed multiplicities and the ratios of the particles in heavy-ion collisions. 



3. Quasiparticle Models (QPM) 



The quasiparticle model in QCD is a phcnomcnological model which is based on the 
idea that quark gluon fluid can be expressed in terms of quasiparticles. It was first 
proposed by Golviznin and Satz [6] and then by Peshier et. al. [7]. In quasiparticle 
models, the system of interacting masslcss quarks and ghions can be effectively described 
as an ideal gas of 'massive' noninteracting quasiparticles. It was assumed that energy w 
and momentum k of the quasiparticles obey a simple dispersion relation : 



u?{k,T) = k'^ +m''{T), 



(2) 



where m(T) is the temperature-dependent mass of the quasiparticle. The pressure and 
energy density of the ideal gas of quasiparticles are dependent on lj and m(T) and are 
given by [4]: 



Pid{T, m) = T 
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where d represents the degeneracy factor for the quarks and/or gluons. However, Goren- 
stein and Yang pointed out that this model involves a thermodynamical inconsistency 
[4]. It did not satisfy the fundamental thermodynamic relation: e(T) = T^^^^ — p{T). 
To remove this inconsistency, two different apporaches were proposed. 

In the first quasiparticle model (QPM I), after reformulating the statistical mechanics 
and incorporating the additional medium contribution, the pressure p and energy density 
e for a system of quasiparticles can be written in a thermodynamically consistent manner 
as follows [4,7]: 

piT,m)=p,d-B{T,^iB), (5) 

€{T,m)=€id + B{T,^iB). (6) 

The first term on the right hand side of both the equations are the standard ideal gas 
expressions given by Eq (3) and (4), respectively. The second term represents the medium 
contribution and is represented by the T and /is dependent bag constant [9]. 

In QPM II, Bannur has suggested that the relation between pressure and grand canon- 
ical partition function cannot hold good if the particles of the system have a temperature- 
dependent mass [5] . So he started with the definition of average energy and average num- 
ber of particles and derived all the thermodynamical quantities from them in a consistent 
manner. In this model, pressure of the system aX jjq = can be obtained as [9] : 

piT,fx, = 0) _ po J eiT,f,,=0) 

To 

where po is the pressure at a reference temperature Tq. We have used Po=0 at To=100 
MeV in our calculation. Using the relation between the number density Ug and the grand 
canonical partition function, we can get the pressure for a system at finite /is [9]. 

4. Results and Discussion 



In order to demonstrate that both types of quasiparticle models reproduce the lattice 
results with the value of parameters chosen here, we find that the variation of pressure 
and energy density with temperature at different Hb/T almost matches with the lattice 
calculations [9] . Additionally, we have plotted here the variation of normalized net quark 
density with respect to temperature at different Hq/T by using QPM II in Fig. 1(a). We 
have used Tq = 100 MeV and At = 115 MeV for a reasonable fit to the data. 

In Fig. 1(b), we have shown the phase boundary obtained in our model. Surprisingly 
we find here that tlic^ first-order deconfining phase transition line ends at a critical end 
point (CEP) and the coordinates of CEP are {Tcep = 183MeV, fxcEP = 166MeV) in 
QPM I and (Tcep = 166MeV,ncEP = 155MeV) in QPM II. It is interesting to find 
that the critical end points obtained by us lie closer to CEP of some lattice calculation 
[16]. We also find a crossover region existing beyond the critical end point where HG 
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Fig. 1. (a) Variation of riq/T''' witii temperature at different Hq/T in Quasiparticlc model QPM 
II. Different symbols are lattice results from Ref. [15]. (b) The location of QCD critical point 
in QCD phase diagram. Pi is the phase boundary in bag model (BM), P2 is the phase bound- 
ary in QPM I and P3 is the phase boundary in QPM II. F is the chemical freezeout line ob- 
tained using our HG model. Ci{Tcep = IGOMeV, hcep = 156MeV) is the CEP on Pi obtained 
in BM, C2{TcEP = 183MeV, hcep = 166Mey) is the CEP on P2 obtained in QPM I and 
CsiTcEP = lQQMeV,ncEP = l5bMeV) is the CEP on P3 obtained in QPM II. 

pressure which is solely dominated by mesonic pressure term in Eq.(l), is always less 
than the QGP pressure. Therefore, no phase transition exists in this region. Since the 
temperature is much higher, the thermal fluctuations break mesonic constituents of HG 
into qtiarks, antiquarks and gluons. 
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